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COLOR  REPRODUCTION  BASED  ON  RED,  GREEN,  AND  BLUE  PRIMARIES 
FOR  A  CYAN-,  MAGENTA-,  AND  YELLOW-BASED  HARDCOPY  DEVICE 

> 

INTRODUCTION 

The  Navy  is  developing  a  database  of  scanned  aeronautical  chart  images,  the  Compressed 
I  Aeronautical  Chart  (CAC),  for  use  in  aircraft  digital  moving-map  systems  and  for  mission  plan¬ 

ning  [1].  The  CAC  uses  a  set  of  30  custom  color  palettes.  Each  palette  consists  oi  240  distinwt 
colors,  and  each  color  is  comprised  of  red-green-blue  (RGB)  intensities.  Intensity  levels  range  from 
0  (no  intensity)  to  2SS  (maximum  intensity).  Although  the  CAC  is  primarily  destined  for  video 
display,  the  need  for  high-resolution  color  hardcopies  also  exists.  Available  plotter  hardware  and 
software  produced  unacceptable  colors,  since  the  low  intensities  of  certain  CAC  palette  colors 
^  resulted  in  poor  color  reproduction  in  hardcopy.  Low-intensity  palette  colors  suffered  from  color  drop¬ 

out,  which  tends  to  produce  colors  with  a  gray-shade  appearance.  It  was  discovered  that  color 
drop-out  is  caused  by  the  addition  of  black  ink,  which  is  specified  by  the  plotter  software.  Some 
color  hardcopy  devices  produce  black  by  blending  the  three  primary  inks  (RGB)  or  cyan,  magenta, 
yellow  (CMY).  However,  this  process  often  results  in  a  somewhat  muddy  black.  This  problem  has 
I  been  solved  for  other  devices  by  using  black  ink  to  produce  true  black.  These  devices  also  add 

black  ink  to  colors  other  than  true  black,  which,  when  combined  with  the  subsequent  reduction  in 
the  CMY  intensities  for  those  colors,  can  produce  the  graying  effect  of  color  drop-out. 

This  report  presents  a  procedure  that  eliminates  color  drop-out  and  reproduces  custom  CAC 
palette  colors  that  are  comparable  in  quality  to  the  colors  found  in  original  aeronautical  charts. 

» 

PROBLEM 

The  custom  color  palettes  for  CAC  data  are  based  on  intensities  of  the  additive  primary  colors;  RGB. 
The  in-house,  high-resolution,  color  plotter  hardware,  which  relies  on  the  complementary  subtractive 
^  primary  colors,  CMY,  produced  plots  with  unacceptable  colors.  In  particular,  the  lower  intensity 

colors  suffered  from  color  drop-out,  which  resulted  in  a  gray-shade  appearance  (Fig.  1).  Higher 
intensity  colors  suffered  less  from  color  drop-out  but  still  did  not  reproduce  well. 

CAC  RGB  intensities  are  normalized  for  use  with  conventional  display  devices  by  using  the 
following  algorithm: 

*  R  =  R/255.0 

G  =  G/255.0 

B  =  B/255.0  (1) 

The  resulting  RGB  values  range  from  0.0  to  1.0  instead  of  0  to  255. 
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Fig.  1  —  CAC  custom  color  palette  with  color  drop-out. 


The  in-house  plotter  plots  CAC  data  by  first  converting  the  normalized  CAC  RGB  intensities  to 
their  CMY  equivalents.  The  following  transformation  [2]  converts  from  RGB  to  CMY: 


c 

1 

R 

M 

= 

1 

- 

G 

Y 

I 

B 

CMY  values  then  are  adjusted  to  reflect  the  plotter’s  addition  of  black  ink.  The  adjustment 
algorithm  [2]  uses  black  (K)  in  place  of  equal  amounts  of  CMY  as  follows: 


K  =  minimum  (C,M,Y) 

C  =C-K 
M  =  M-K 

Y  =  Y  -  K  (3) 


Using  a  4  X  4  pixel  pattern,  CMYK  inks  are  deposited  on  paper  as  a  grid  of  colored  dots.  The 
orientation  of  each  colored  dot  allows  the  eye  to  spatially  integrate  light  that  is  reflected  from 
adjacent  dots.  For  any  given  CMYK  component,  a  7%  minimum  intensity  is  required  to  have  at 
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least  1  element  activated,  out  of  16  in  the  pixel  matrix.  CMYK  components  with  intensities  of  less 
than  7%  (i.e.,  no  color)  will  not  be  plotted.  The  following  algorithm  [3]  is  used  to  calculate  the 
number  of  elements  (within  the  matrix)  to  activate  for  a  given  CMYK  intensity; 

Number  of  elements  =  (intensity  *  matrix  size)  /  100  ,  (4) 

where  intensity=  (intensity  x  100) +  0.5,  and  matrix  size  =  4x4  elements. 

The  CAC  palette  has  been  sorted  by  increasing  intensity.  As  shown  in  Fig.  1,  color  drop-out 
was  most  severe  in  the  lower  intensities.  One  attempt  to  alleviate  color  drop-out  was  based  on  the 
hue,  lightness,  saturation  color  model  [2]  in  which  the  CAC  color  palette  was  shifted  toward  white 
along  the  achromatic  axis  to  increase  the  overall  color  palette  intensities  (Fig.  2).  This  achromatic 
shift  succeeded  in  producing  less  color  drop-out  in  lower  intensity  colors.  However,  the  shift 
resulted  in  the  opposite  problem:  lighter  intensity  colors  were  “washed  out"  (i.e.,  too  much  white). 


Myrick,  Lohrenz,  and  Wischow 


Applying  a  combination  of  achromatic  shifts  to  individual  colors  was  deemed  undesirable  due  to 
the  significant  number  of  colors  that  would  have  to  be  modified  (240  colors  in  each  of  the  30 
standard  CAC  color  palettes).  Even  with  the  availability  of  examples  of  computer-generated 
colors  [4],  where  percentages  of  each  CMYK  ink  are  provided,  there  would  be  simply  too  many 
individual  CAC  palette  colors  to  be  matched. 

Another  attempt  to  alleviate  color  drop-out  involved  increasing  the  pixel  pattern  to  8  x  8.  This 
procedure  required  a  minimum  intensity  of  2%  for  any  given  CMYK  color  to  have  at  least 
1  element  activated  out  of  64.  Although  the  resulting  plots  had  less  color  drop-out,  this  attempt 
proved  to  be  unsatisfactory  because  they  suffered  a  gridded  appearance  due  to  the  orientation  of 
colored  dots  within  the  matrix. 

An  underlying  problem  with  both  of  these  approaches  is  the  large  number  of  colors  to  be 
manipulated.  Rather  than  working  with  every  color  in  a  given  palette,  a  test  palette  was  devised. 
Since  CAC  data  is  comprised  of  RGB  intensities  (which  are  later  transformed  into  their  CMYK 
components),  various  intensities  of  RGB  were  selected  as  test  colors.  The  intensities  of  the  colors 
were  chosen  to  vary  in  lightness  by  12%,  from  little  to  full  pure  color.  Figure  3  presents  a  plot  of 
these  colors  and  shows  how  adding  black  ink  influences  the  appearance  of  the  resultant  RGB 
colors.  Except  in  the  case  of  colors  with  full  intensity  (e.g.,  pure  red,  pure  green,  or  pure  blue), 
black  ink  is  always  noticeably  present. 


SOLUTION 

The  CAC  color  palettes  contain  true  black,  but  true  black  is  rarely  found  in  the  digitized 
aeronautical  charts.  It  is  present  in  CAC  palettes  primarily  for  areas  of  no  data  coverage.  Since  the 
CAC  data  rarely  (if  ever)  require  true  black,  and  since  the  addition  of  black  by  the  plotter  seriously 
compromises  the  quality  of  many  other  CAC  colors,  the  final  solution  to  the  color  drop-out  problem 
pivoted  around  the  elimination  of  black  ink  from  the  plotting  process. 

The  FORTRAN  source  code  for  the  plotter  software  [3]  was  modified  in-house  to  completely 
omit  the  inclusion  of  black  ink:  RGB  intensities  were  converted  to  their  CMY  equivalents  (Eq.  2), 
but  the  adjustment  algorithm  (Eq.  3)  for  black  ink  was  eliminated.  Figure  4  presents  the  resulting 
test  plot  of  RGB  colors  (in  which  intensity  varies  by  12%  from  little  to  pure  color)  using  CMY  inks 
and  no  black  ink.  The  lower  intensity  RGB  colors,  which  previously  suffered  from  color  drop-out, 
now  are  acceptable  in  appearance.  The  higher  intensity  RGB  colors  are  also  acceptable.  Using  the 
modified  plotter  software,  another  custom  color  palette  plot  was  created.  As  shown  in  Fig.  5,  the 
palette  suffered  no  appreciable  color  drop-out,  and  all  colors  were  more  accurately  reproduced. 


CONCLUSIONS  AND 
RECOMMENDATIONS 

Results  from  the  color  experiments  documented  here  indicate  that  the  addition  of  black  ink 
interferes  with  the  reproduction  of  custom  colors  on  some  hardcopy  plotters.  The  Navy  standard 
CAC  database  was  used  in  this  study.  CAC  data  use  a  standard  set  of  30  custom-designed  color 
palettes,  each  of  which  consists  of  240  RGB  colors.  The  custom  RGB  colors  in  CAC  data  did  not 
reproduce  well  on  an  in-house,  CMYK-based,  color  hardcopy  device  that  added  black  ink  during 
the  RGB  to  CMYK  transformation.  In  particular,  the  addition  of  black  ink  resulted  in  a  graying 
effect,  termed  color  drop-out,  in  most  CAC  colors.  Several  adjustments  were  considered,  including 
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Fig.  S  —  CAC  custom  color  palette  without  color  drop-out. 


Fig.  6  —  Reference  plot  with  CAC  data. 
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a  shift  of  the  colors  along  the  achromatic  axis  (toward  higher  overall  intensities),  increasing  the  size 
of  the  pixel  pattern  matrix  from  4  x4  to  8  x  g  pixels,  and  eliminating  the  inclusion  of  black  ink. 
The  first  two  approaches  did  not  adequately  solve  the  color  reproduction  problems;  the  third  approach 
proved  to  be  the  best  solution.  Since  true  black  is  rarely  needed  by  the  CAC  database,  and  since 
black  can  be  approximated  by  blending  CMY  colors,  the  use  of  black  ink  was  completely  eliminated. 
Test  plots  revealed  that  the  deletion  of  black  ink  from  the  RGB  to  CMYK  transformation  resulted 
in  the  successful  reproduction  of  all  CAC  custom  colors.  Appendix  A  contains  source  code  for  the 
new  program,  DISPLAY_CAC_CALCOMP.FOR.  that  is  used  to  plot  CAC  data  on  the  electrostatic 
plotter. 

An  interesting  corollary  to  this  study  was  the  demonstrated  vulnerability  of  plot  quality  to 
environmental  and  chemical  factors.  The  quality  of  plot  appearance  was  profoundly  influenced  by 
environmental  factors,  such  as  room  temperature  and  humidity.  Chemical  factors  relating  to  the 
toner  and  replenisher  (i.e.,  age  and  percentage  used)  also  affected  plot  quality.  The  adverse  influ¬ 
ence  of  these  factors  was  manifested  by  uneven  distribution  and  absorption  of  ink.  Due  to  the 
significant  influence  of  environmental  and  chemical  factors  on  plot  quality  and  appearance,  reference 
plots  are  now  included  with  ail  in-house  plotter  output.  These  reference  plots,  which  include  smaller 
versions  of  the  RGB  and  color  palette  plots,  help  to  identify  whether  a  poor  plot  was  due  to  color 
problems  or  to  environmental/chemical  influences.  A  typical  reference  plot,  along  with  a  portion 
of  CAC  data,  is  shown  in  Fig.  6.  Appendix  B  lists  source  code  for  the  subroutine, 
REFERENCE_PLOT.FOR,  that  is  used  to  generate  reference  plots. 
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APPENDIX  A 


TITLE:  Di splay _CAC_CalComp. FOR 


DESCRIPTION; 


This  program  creates  a  plot  file  of  Compressed 
Aeronautical  Chart  (CAC)  data.  Specified  rows  and 
columns  of  CAC  segments  are  decompressed  and  written 
to  a  CalComp- formatted  file.  This  file  contains  the  required 
format  for  plotting  CAC  segments  on  a  CalComp  electrostatic 
color  plotter. 

The  CAC  is  produced  by  the  Naval  Research  Laboratory  (NRL) , 
Stennis  Space  Center,  Mississippi.  Some  familiarity 
with  CAC  is  required  to  execute  this  program. 

See  NOARL  Report  8,  July  1990  (product  specification) 
for  additional  information  about  CAC. 

Reference  plots  which  include  a  Red,  Green,  Blue  color  scale 
and  a  CAC  color  palette  are  optional. 

OUTPUT:  A  binary  plot  file  (in  CalComp  format)  containing  CAC 

coverage  over  the  specified  rows  and  columns  is  written  in  the 
user's  current  directory.  The  file  name  is  supplied  as 
program  input . 

AUTHOR:  Stephanie  My rick 
NRL  Code  7441 

Stennis  Space  Center,  MS  39529-5004 


program  DISPLAY_CAC_CALCOMP 
implicit  none 

include  ' ($RMSDEF) '  !  VAX/VMS  System  Definition 

byte  rgbtable(256, 3)  !  Color  table  extracted  from  color 

!  "palette. dat" 

byte  seg_Duf fer (256, 256)  1  Uncompressed  segment 

byte  iktmp(4)  !  Temp  variable  equivalenced  with  KITMP 
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character*80 

character *  80 

character* 132 

character*132 

character* 132 

character*80 

character*80 

character*! 

character*4 

character*6 

character*6 

character*! 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

logical*! 

logical*! 

real*4 

real*4 

real*4 

real*4 

real*4 

real*4 


plotf ilename 

filename 

pal name 

wild_card 

pa_f ilename 

pa_dir 

pa_f ound 

scale_char 

pa_char 

row_char 

col_char 

reply 


Output  CalComp  plot  file  name 
Compressed  segment  file  name 
Full  color  palette  file  name 
Wild  card  file  specification 
PA  directory  as  a  file  name 
Palette  directory  path 
Palette  number  found 
Chart  scale 
Palette  number 
Row  number,  user  input 
Column  number,  user  input 
User  reply,  input 


length  !  Length  of  plotfilename 

seg_plot (256 , 256)  !  Compressed  segment  plot  buffer 


kitmp 

row 

col 

start_row 
start_col 
end_row 
num_rows 
num_cols 
markl , mark2 
pa_dir_len 
pa_f ound_len 
pal_color 
palette_of f set 
pal_shif t 
Idev 

f ile_stat 

f ile_len 

context , status 

input_stat 

isize 

zone 


Temp  variable  equiv.  with  IKTMP 
Current  row 
Current  column 
Start  row  of  display 
Start  col  of  display 
Last  row  of  display 
Number  of  row  to  display 
Number  of  columns  to  display 
Denotes  substring  locations 
Length  of  PA  directory  name 
Length  of  palette  name 
Color  palette  index 

Offset  for  loading  our  color  palette 
Palette  shifting 
Logical  device  number 
File  open  status 

Length  of  compressed  segment  file 
Used  w/RMS  calls  to  locate  files 
Status  of  user  input 
COLRGN  argument,  matrix  size 
M4  zone 


i,j  !  Temporary  array  indeces 

DECOMPRESS_SEGMENT  !  Routine  to  open  &  read 

!  a  compressed  segment 

LIB$FIND_FILE  !  RMS  routine  to  locate  a  file 


want_ref 

NO_MORE_FILES 

red, grn, blu 

x_axis,y_axis 

x_position 

y_position 

ifactor 

ifactor  offset 


Indicates  use  of  RGB  and 
color  palette  reference  plots 
RMS  file  search 

RED  GREEN  BLUE  colors 

X-  and  y-axis  plot  origins 

x-axis  plotting  position 

y-axis  plotting  position 

plot  expansion  factor 

factor  dependent  offset  for  plots 


equivalence  (kitmp, iktmp) 
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C*  Initializations  * 

if  it  it  it  it  it  it  "k  it  it  if  it  if  ^  'it  it  it  it  it  'it 


Idev  =  9  ! 
isize  =  4  ! 
x_axis  =0.0  ! 
y_axis  =0.0  ! 
x_position  =0.0  ! 
y_position  =0.0  ! 


Output  file  logical  device  number 
Use  this  as  the  default  size 
X  axis  origin 
y  axis  origin 

initial  x-axis  plot  position 
initial  y-axis  plot  position 


call  INIT_MEM  !  Initialize  memory  for  CAC  utility  software 


C****  Open  output  plot  file  **** 
write  ( 6 , * )  '  ' 

write  (6,*)  'Enter  the  output  CalComp  plot  file  name' 

write  (6,*)  '{Omit  a  file  extension)' 

write  (6,*)  '  ' 

read  (5, '(a)')  plotf  ileneime 

call  STRING_LENGTH  (plotfilename, length) 

plotfilename(length+l:length+4)  =  '.DAT'  I  the  default  extension 

open  (Idev, f ile=plotf ilename, status='new' , 

&  recordtype= 'variable' ,  blocksize=484 , err=9001 ) 

write  (6,*)  '  ' 

write(6,*)  'Opened  CalComp  plot  file:  ', plotfilename (1 : length+4 ) 


C****  Calcomp  call:  plot  initialization  **** 
call  plots (x_axis,y_axis, Idev) 


C****  returns  here  via  ''Z  to  re-enter  program  input  **** 

100  write  (6,*) 

input_stat=-l 

do  while  (input_stat.  ne.  0) 

write(6,'("  Enter  SCALE  (CTRL_Z  to  quit):  '',$)') 
read(5, ' (a) ' , iostat=input_stat, end=9000)  scale_char 
enddo 


write  (6,*) 
input_stat=-l 

do  while  (input_stat  .ne.  0) 

write(6,'("  Enter  PA#:  ",$)') 
read(5, ' (a) ' , iostat=input_stat , end=9000) 
enddo 


pa_char 


write  (6,*) 
input_stat=-l 

do  while  (input_stat  .ne.  0) 

write(6,'(''  Enter  ZONE  (0,1, 2, 3, 4):  '',$)') 
read ( 5 , * , iostat=input_stat , end=9000 )  zone 
enddo 


If’ 
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write  (6,*)  '  ' 
input_stat=-l 

do  while  (input_stat  .ne.  0) 

write(6, *) 'Do  you  want  to  include  the  RGB  and' 
write(6,'(''  color  palette  reference  plots?  (Y/N)  '',$)') 
read ( 5 , ' (a) ' , iostat=input_stat , end=9000 )  reply 
enddo 

C****  By  default,  the  reference  plots  are  included  **** 
if  ((reply  .eq.  'Y')  .or.  (reply  .eq.  'y')  .or. 

St  (reply  .eq.  '  ')  )  want_ref  =  .true. 

C****  you  need  to  set  IFACTOR  depending  on  the  isize  selected  **** 

C****  plotting  offset  will  not  change  **** 

ifactor  =  0.5 
ifactor_of f set  =  5.12 

C****  construct  directory  path  **** 

pa_dir  =  ' CHART_ODI_DISK: [MAP'  //  scale_char  //  ']' 
call  STRING_LENGTH  (pa_dir , pa_dir_len) 

pa_ found  =  '  ' 
pa_f  ilenaitie  =  '  ' 
context  =  0 
wild_card  =  '*.*;*' 

NO_MORE_FILES  =  .false. 

C****  search  for  palette  directories  **** 

do  while  ( (pa_char  .ne.  pa_found(4 : 6) )  !  palette  #'s  don't  match 

St  .or.  (NO_MORE_FILES  .eq.  .false.)) 

status  =  LIB$FIND_FILE( (  pa_dir ( 1 : pa_dir_len)  // 

'pa* .dir; * ' ) , 

pa_f ilename,  context,  wild_card) 
if  (status  .eq.  RMS$_NMF)  then 
continue 

else  if  (status  .eq.  RMS$_NMF)  then 
NO_MORE_FILES  =  .true. 

else 

NO_MORE_FILES  =  .true, 
end  if 

MARKl  =  index  (  pa_f ilename,  ']') 

MARK2  =  index  (  pa_f ilename,  '.DIR') 
pa_found  =  pa_f ilename (MARKl +1 :MARK2-1 ) 


enddo 


(?’  Ri  R>  R> 
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C****  Read  the  color  palette  and  separate  the  RGB  values  **** 

•  C****  for  loading  the  color  table  into  the  graphics  processor  **** 

call  STRING_LENGTH  (pa_found, pa_f ound_len) 

pa_dir  =  pa_dir ( 1 :pa_dir_len  -  1)//  // 

pa_found ( 1 :pa_found_len)  // 

'  ]  ' 

call  STRING_LENGTH  (pa_dir,  pa_dir_len) 

READ_PALETTE  is  a  C  function,  hence  character  string  **** 
must  be  terminated  with  a  null  character  **** 

palnaime  =  {pa_dir  ( 1  :pa_dir_len) //' palette . dat '  //  char(O)  ) 


call  READ_PALETTE  ( %val ( 0 ) , %ref (palname ) , 

rgbtabled,  1)  , 
rgbtable (1, 2) , 
rgbtable (1,3) ) 


C 

C 

C 

C 


it  ir  ie  it 
it  ir  it  it 
it  ir  it  "k 
it  it  it  it 


Load  the  color  palette:  RED,  GREEN,  BLUE 

An  offset  of  250  is  required  to  skip  over  memory  locations 
that  are  reserved  for  CalComp  pre-defined  symbols  in  use 
when  isize  =4.  We  will  use  an  offset  of  255 
palette_of f set  =  255  !  255  for  isize  =  4 

!  0  for  isize  =  8 

do  i  =  1,240 


*  *  *  * 
it  it  it  it 

*  *  it  it 

*  it  *  it 


IKTMP(l)  =  rgbtabled,  1) 
red  =  kitmp 
red  =  red/255.0 


!  used  to  be  +  pal_shift 
!  remove  declaration  if  OK 


iktmpd)  =  rgbtable  ( i ,  2 ) 
grn  =  kitmp 
grn  =  grn/255,0 

iktmpd)  =  rgbtable (i,  3) 
blu  =  kitmp 
blu  =  blu/255.0 

C  ****  Calcomp  call:  generate  a  user-defined  color  **** 

call  colrgn(2, red, grn, blu, i+palette_of f set , isize) 


enddo 


C****  Determine  ROW/COL  values  of  data  in  specified  directory  **** 
Qit***  Return  here  via  ''Z  to  re-enter  ROW/COL  info  *★** 

150  write  (6,*) 

input_stat=-l 

do  while  (input_stat  .ne.  0) 

write(6,'(''  Enter  start  ROW  (CTRL_Z  to  change 


o  o 
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&  ROW/ COL) : 

read(*, ' (a) ' , iostat=input_stat, end=150)  row_char 
i f  ( row_char  . ne .  '  ' )  then 

read  (row_char, ' (bn, i6) ' )  start_row 
end  if 
enddo 

write  ( 6 , * ) 
input_stat=-l 

do  while  (input_stat  .ne.  0) 

write(6,'{''  Enter  start  COL  (CTRL_Z  to  change 
&  ROW/COL)  :  ",$)') 

read(*, ' (a) ' , iostat=input_stat, end=150)  col_char 
if  (col_char  .ne.  '  ')  then 

read  (col_char, ' (bn, i6) ' )  start_col 
end  if 
enddo 

write  (6,*) 
input_stat=-l 

do  while  (input_stat  .ne.  0) 

write(6,'(''  Enter  the  number  of  ROWs  to  plot 
&  (CTRL_Z  to  change  ROW/COL): 

read(*, ' (a) ' , iostat=input_stat, end=150)  row_char 
i f  ( row_char  . ne .  '  ' )  then 

read  (row_char,  '  (bn,  i6) ' )  n\im_rows 
end  if 
enddo 

write  (6, *) 
input_stat=-l 

do  while  (input_stat  .ne.  0) 

write(6,'(''  Enter  the  number  of  COLs  to  plot 
&  (CTRL_Z  to  change  ROW/COL) : 

read(*, ' (a) ' , iostat=input_stat, end=150)  col_char 
if  (col_char  .ne.  '  ')  then 

read  (col_char, ' (bn, i6) ' )  num_cols 
end  if 
enddo 

call  draw_infobox  (palncime,  scale_char,  zone,  start_row, 

*  start_col,  num_rows,  num_cols) 

C****  If  desired,  display  RGB  and  color  palette  reference  plots  **** 
if  (want_ref)  then 

call  reference_plot 

x_position  =5.0  !  Reference  plot  has  used  approx.  5  inches, 

C  *****  calComp  call:  adjust  X  and  Y  origins  ***** 

call  plot  (x_position,y_position, -3 ) 
endif  !  adjust  x-axis 


n  n  o  n  n  n  n 


j 
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C 


•k  ie  it  if  * 


CalComp  call:  expand  plot  by  this  factor 
call  factor ( ifactor) 


it  it  it  it  k 


Q*-k***,  Need  to  replace  the  directory  terminator  ']'  with  ***** 
C*****  a  subdirectory  ' . '  in  PA_DIR  for  appending  Row  ***** 
C*****  subdirecotry  specification  .Rsnnnnn]  ***** 

•  call  STRING_LENGTH (pa_dir, pa_dir_len) 

pa_dir  =  pa_dir(l:pa_dir_len-l)  //  //  char(O) 


C****  Display  compressed  chart  data  plots  with  the  upper  right  segment  **** 


end_row  =  start_row  +  num_rows  -1 

do  row  =  end_row, start_row, -1  !  each  iteration  plots  a  row 

do  col  =  start_col , start_col+num_cols-l  !  each  iteration  plots  a  col 

*****  Build  file  name:  ***** 

*****  get_SEGMENT_NAME  is  a  C  function.  Hence  character  ***** 
*****  strings  must  be  terminated  with  a  null  character  ***** 
*****  arguments  must  be  passed  by  value  &  reference  ***** 

call  GET_SEGMENT_NAME (%ref (pa_dir) , %val (row) , 

Sc  %val  (col)  ,  %val  (zone)  , 

&  %ref (filename)  ) 


C 


7000 


call  STRING_LENGTH  ( filename, file_len) 

****  Read  and  decompress  a  segment  file.  If  no  **** 

****  file  exists,  a  "blank"  segment  (i.e‘.  a  **** 

****  hole)  in  the  plot  results  **** 

file_stat  =  DECOMPRESS_SEGMENT  ( 

%ref ( filename) , 
seg_buf fer) 

****  Check  for  irregular  return  status,  normal  =  1  **** 
if  {file_stat  .ne.  1)  then 
print* 

print*, '»  No  data  for  ROW', row,  '  COL  ',col 
print*,'  ' ,  f  ileneune  ( 1 :  f  ile_len)  ,  '  NOT  found...' 
else 

write (6, * ) 

write(6,7000)  row,  col 

f ormat ( '  Loading  data  for  ROW:  ',i6.5,'  COL:  ',i6.5) 


ff’ 
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C 


C 


*****  load  data  for  plotting  ***** 
do  i=l,256 
do  j=l,256 

IKTMP(l)  =  seg_buffer(i, j) 

seg_plot ( i , j )  =  KITMP  +  palette_of f set  +  1 

*****  Check  for  colors  out  of  range  ***** 
if  ( (kitmp+palette_of f set+1  .It.  256)  .or. 
(Kitmp+palette_of f set+1  .gt.  495)) 
write (6,*)  'Indexing  color  out  of  range', 
kitmp+palette_off set+1 


enddo 

enddo 

endif  !  for  file_stat  =  1 


C 


*****  Calcomp  call:  adjust  x  &  y  plot  positions 

call  plot  (x__position,y_position,  3 ) 

type  *, 'plotting  x  Se  y ' , x_position,y_position 


C 


*****  CalComp  call:  fill  plot  buffer  ***** 
call  rasfil (0.02, 0 . 02, 256, 256, seg_plot) 


C  *+**  adjust  y  coordinate  -  move  over  250  pixels  **** 

C  *  +  *★  for  plotting  the  next  segrment/ column  **** 

y_position  =  y_position  +  ifactor_of f set 

enddo  !  EndDo  columns 


C  ****  Adjust  x  coordinate  -  move  up  250  pixels  **** 

C  ****  for  plotting  the  next  row  of  segments  **** 

x_position  =  x_position  +  ifactor_of f set 

C  **★*  Adjust  y  coordinate  -  reset  to  origin  **** 

C  **★*  for  plotting  the  next  column  of  segments  **** 

y_position  =  0.0 

enddo  !  EndDo  rows 

C  *****  Calcomp  call:  fill  every  image  pixel  ***** 

C  call  rasfil  (C  ,01, 1000, 1000,  seg_plot) 


C 


★  ★  ★  ★  ♦ 


CalComp  call:  enlarge  the  size  of  the  entire  plot 


it  it  -k  *  It 


c 


9000 
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*****  Calcomp  call:  terminate  plot  ***** 
call  plot ( 0 0 999 ) 

close ( Idev) 

stop  '  End  of  program  DISPLAY_CAC_CALCOMP' 

stop  '  Error  opening  CalComp  plot  file!  Aborting  the  program' 


9001 


APPENDIX  B 


c 

C  TITLE:  Ref erence_Plot . FOR 
C 

C  DESCRIPTION:  A  FORTRAN  subroutine  which  utilizes  CalComp  software  for 
C  plotting  Red,  Green  and  Blue  (RGB)  intensities  and  a 

C  CAC  custom  color  palette. 

C  20  RGB  Intensities,  varying  by  12.75%  simulate  the  original 

C  CAC  intensities  that  range  from  0-255. 

C 

C  ARGUMENTS:  None 
C 

C  RETURNS :  None 
C 

C  REQUIRED  SUBROUTINES: 

C 

LOAD_PLOT_BUF :  Fills  the  plot  buffer  with  25x25  pixel  "blocks" 

of  color.  The  array  dimensions  of  the  plot  buffer 
are  adjustable. 


C 
C 
C 
C 
c 
c 
c 
c 

subroutine  reference_plot 


AUTHOR:  Stephanie  Myrick 
NRL,  Code  351 

Stennis  Space  Center, MS  39529-5004 


implicit  none 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 

integer*4 


rgb_buf (500,500)  !  Plot  buffer 
palette_buf (400, 400)  !  Plot  buffer 


pal_color 
num_horiz 
num_vert 
num _pixels 

rgb_of f set 
palette_of f set 

isize 

i 


!  Color  palette  index 
!  Number  of  horizontal  color  blocks 
!  Number  of  vertical  color  blocks 
!  Number  of  pixels  comprising  the 
!  color  block  dimensions,  x-y  plane 
!  Offset  for  loading  RGB  values 
!  Offset  for  loading  our  color 
!  palette 

!  COLRGN  argument,  matrix  size 
!  DO  loop  counter 


17 


n  o 


Myrick,  Lohrenz,  and  Wischow 


18 


integer*4 

integer*4 

real*4 

real*4 

real*4 

real*4 

real*4 

C*  Begin  * 


x_buf 

y_buf 

x_position 
y_position 
red, grn, blu 
rgb 

intensity 


isize  =  4 


Plot  buffer  x  dimension 
Plot  buffer  y  dimension 

Plot  x-axis  position 
Plot  y-axis  position 
individual  RED  GREEN  BLUE  colors 
RGB  intensity 

incremental  intensity  for  RGB 


C****  Create  pure  RGB  colors  that  vary  in  20  intensities  **** 
C****  Intensities  range  from  little  to  pure  color  **** 
C****  Allow  enough  offset  for  reserved  locations  &  palette  **** 
C****  where  250  locations  are  reserved  for  CalComp  symbols  **** 
C****  and  240  locations  are  reserved  for  CAC  color  palette  **** 


rgb  =  12.75 

rgb_offset  =  500  !  memory  location  offset 

intensity  =  12.75 

do  i  =  1,19 

red  =  rgb  /  255.0 
grn  =  rgb  /  255.0 
blu  =  rgb  /  255.0 

call  colrgn  (2 , red, 0 . 0, 0 . 0, i+rgb_of f set , isize) 
call  colrgn  (2, 0.0,grn,0.0, i+rgb_offset+20, isize) 
call  colrgn  (2 , 0 . 0 , 0 . 0,blu, i+rgb_of f set+40 , isize) 

C  **★*  Increase  color  intensity  by  l/20th  **** 

rgb  =  rgb  +  intensity 


enddo 


*****  load  pure  RGB  values  ***** 

*****  since  255  =  pure  color,  pass  1  as  intensity  ***** 
call  colrgn  (2 , 1 . 0 , 0 . 0 , 0 . 0 , 20+rgb_of f set , isize) 
call  colrgn  (2 , 0 . 0, 1 . 0, 0 . 0, 40+rgb_of fset , isize) 
call  colrgn  (2 , 0 . 0 , 0 . 0, 1 . 0 , 60+rgb_of fset , isize) 

Q****  Fill  a  3x20  color  matrix  with  the  RGB  intensities  ***** 
num_horiz  =  3  !  Each  horizontal  is  an  RGB  primary 

num_vert  =  20  !  These  contain  the  intensities 
num_pixels  =  25  !  Each  block  of  color  is  25x25  pixels  wide 


x_buf  =  500  !  RGB_buf  array  x-dimension 


(fi  Ri  R’  R’ 
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» 


i 


y_buf  =  500  !  y-dimension 

call  load_ploc_buf  (nuin_horiz, nuin_vert , 

nuin_pixels,  rgb_of  f  set ,  rgb_buf , 
x_buf ,y_buf ) 

Initialize  x-axis  and  y-axis  positions  **** 
x_position  =  0.0 
y_position  =  0.0 

call  plot (x_position,y_position, 3) 

CalComp  call  to  fill  the  raster  image  **** 
call  rasfil  (0 . 01, 0 . 01, 500, 500, rgb_buf ) 

CalComp  call  to  modify  plotting  coordinates  **** 

Need  to  move  over  for  the  next  plot  **** 

x_position  =  0.0 
y_position  =  2.0 

call  plot  {x__position,y_position,  3) 

Load  the  CAC  color  palette.  An  offset  of  250  is  required  **** 
to  skip  over  memory  locations  that  are  reserved  for  CalComp  **** 
pre-defined  symbols  in  use  when  isize  =4.  we  will  use  255  **** 
palette_of f set  =  255 

Fill  a  15x16  matrix  of  240  CAC  palette  colors  **** 
num_horiz  =  15 
num_vert  =  16 
x_buf  =  400 

y_buf  =  400  , 

call  load_plot_buf  (num_horiz, n\im_vert , 

num_pixels, palette_of f set , palette_buf , 
x_buf , y_buf ) 

C****  CalComp  call  to  fill  the  raster  image  **** 

call  rasfil  (0 . 01, 0 . 01, 400, 400,palette_buf ) 

return 

end 


I 
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C 

C  TITLE:  LOAD_PLOT_BUF 
C 

C  DESCRIPTION:  Fills  the  buffer  containing  reference  plots. 

C  Creates  a  matrix  of  square  color  "blocks"  which  are 

C  composed  of  (NUM_PIXELS  x  NUM_PIXELS)  pixels. 

C  Colors  beginning  at  memory  location  OFFSET 

C  are  used  to  fill  these  "blocks"  of  color. 

C 

C  ARGUMENTS : 

C  num_horiz 

C  num_vert 

C  num_pixels 

C  offset 

C 

C  plot_buf 

C  x-buf 

C  y-buf 

C 

C  RETURNS:  None 

C 

C  AUTHOR:  S. Myrick 


Total  number  of  horizontal  color  "blocks" 

Total  number  of  vertical  color  "blocks" 

Color  "block"  dimension.  #  of  pixels  in  x-y  plane 
Memory  offset,  pointing  to  first 
location  storing  color 

Buffer  containing  reference  plots,  ajustable  array 
Adjustable  array  x  dimension 
Adjustable  array  y  dimension 


subroutine  load_plot_buf  (num_horiz,num_vert,num_pixels, 

&  of fset,plot_buf ,  x_buf,y_buf) 


Total  number  of  horizontal  color  "blocks" 

Total  number  of  veritcal  color  "blocks" 
Dimension  of  color  "block"  in  x-y  planes 
Memory  offset,  pointing  to  first  location 
!  storing  color 

plot  buf  array  dimension 
plot  buf  array  dimension 
.buf,  y_buf)  !  Buffer  containing  reference  plots 
Counter  of  horizontal  color  blocks 
Counter  of  vertical  color  blocks 
Current  pixel  row 
Current  pixel  colijmn 

Starting  row  for  creating  color  "block" 

Starting  column  for  creating  color  "block" 

Last  row  for  creating  color  "block" 

Last  column  for  creating  color  "block" 

Color  used  in  filling  color  "block" 


implicit  none 

integer*4 

num_horiz  ! 

integer*4 

num_vert  ! 

integer*4 

num_pixels ! 

integer*4 

offset  ! 

integer*4 

x_buf  I 

integer*4 

y_buf  ! 

integer*4 

plot_buf  (x_l 

integer*4 

horiz  ! 

integer*4 

vert  ! 

integer*4 

row  ! 

integer*4 

col  i 

integer*4 

startrow  ! 

integer*4 

startcol  ! 

integer*4 

endrow  ! 

integer*4 

endcol  ! 

integer*4 

color  1 

u  u  u 


» 
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C*  Begin  * 

color  =  offset  !  -  1  adjust  for  DO  loop  incrementation 
startrow  =  1 
endrow  =  num_pixels 
startcol  =  1 
endcol  =  nuin_pixels 

C****  Create  a  pixel  image,  or  “block".  **** 

C****  Each  iteration  of  the  outer  loop  creates  NUM_HORIZ  horizontal  **** 
C****  blocks  of  color.  **** 

C****  Each  iteration  of  the  inner  loop  creates  NUM_VERT  blocks  of  color.* 
C****  Hence,  NUM_HORIZ  x  NUM_VERT  blocks  of  individual  color  are  created* 

do  horiz  =  l,num_horiz 


***  Create  “blocks"  of  color.  Each  iteration  fills  *** 
***  one  block  that  is  NUM_PIXELS  pixels  wide  with  one  color  *** 
***  These  blocks  are  built  or  placed  horizontally.  *** 


do  vert  =  1 , num_vert 

color  =  color  +  1 
do  col  =  startcol , endcol 


get  a  color 

I  for  colums/pixels 


do  row  =  startrow, endrow 
plot_buf (row, col)  =  color 
enddo 


!  fill  rows/pixels 
!  with  this  color 


enddo 


startcol 

endcol 


startcol  +  num_pixels  !  move  over  to  fill 
endcol  +  nuin_pixels  !  another  25  columns  of 

!  pixels 


enddo  !  for  tilling  in  a  vertical  block  of  color 


startcol  =  1  !  restart  another  set  of  columns/pixels 

endcol  =  num_pixels 

startrow  =  startrow  +  num_pixels  !  move  over  to  fill 

endrow  =  endrow  +  num_pixels  !  another  row  of  pixels 

enddo  !  for  filling  in  a  horizontal  block  of  color 


return 

end 
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C  + 

it  it  it  -k  -k  it  it  *  it  it  it  it  it  it  it  it  -k  it  ir  -k  it  ir  ii  ir  -k  ii  it  it  it  it  it  it  -k  it  -k  it  it  ii  ir  it  -k  it  it  -k  it  it  it  it  ii  it  -k  it  *  *  *  it  it  it  it  *  lr  it  it  *  -k  *  *  *  * 

C  NAME:  draw_infobox 
C 

C  PURPOSE:  To  draw  an  information  box  using  calcomp  calls  that  contains  the 
C  palette  used,  the  start  row,  the  start  col,  the  number  of  rows  and 

C  cols,  the  scale,  and  the  zone. 

C 

C  AUTHOR:  Joyce  Michelle  Mehaffey  (Planning  System  Inc.),  Nov  92 
C 

C  SUMMARY: 

C  parameter  list  NONE 
C 

C  DESCRIPTION:  This  subroutine 
C  this  routine  is 

C 

c 

C  SYSTEM  DEPENDENT  FUNCTIONS: 

C  When  linlcing  this  program, 

C  =  = 

Qlt^*****-*****«*********************1^***************************it**lk^******* 


assumes  that  a  plot  file  has  been  open  before 
called. 


you  must  link  in  the  calcomp  library 


subroutine  draw_infobox  (pa_dir,  scale_char,  zone,  srow, 

scol,  nrows,  ncols) 


I 


I 


» 


implicit  none 


include  'v2_dir:m4_constants. inc' 


integer*4 

integer*4 

integer*4 

real*8 


scale, zone, slen 
srow,  scol 
nrows,  ncols 

y_pos,  x_pos  !  start  position  for  a  plot 


character*80 

character*80 

character*3 

character*3 

character*! 

character*2 

character*23 


stringl  !  filename  for  the  plot  file 

pa_dir 

srow_char,  nrow_char 
scol_char,  ncol_char 
scale_char 
zone_char 
create_date 


x_pos  =  0.0 
y_pos  =  9.0 

call  plot(x_pos,  y_pos,  3) 

factor  size  =  1.0 
boldface  lines 


call  factor  (1.0) 
call  newpen(ll) 


!  always  make  the 
!  set  pen  to  draw 
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call  plot  (  0.0,16.0,2)  .'  draw  thicJc  outside  border 

call  plot (  1.85,16.0,2) 

call  plot (  1.85,  9.0,2) 

call  plot (  0.0,  9.0,2) 

call  plot (  1.0, 9.1  ,3) 

call  newpen(3)  1  set  pen  to  draw  boldface  lines 

call  lib$date_time  (create_date)  !  creation  date 
call  symbol (1.50,9.25, .10, ' created: ',,90.0,8) 
slen  =  INDEX  (create_date, ' . ' )  -1 

call  symbol (1.50,10.25, .10, create_date, , 90 . 0 , slen) 

call  strSupcase (pa_dir , pa_dir) 

call  STRING_LENGTH(pa_dir, slen) 

call  symbol (0.25,9.25, .12, pa_dir , , 90 . 0 , slen) 

call  symbol ( 0 . 50, 9 . 25 , . 12 , ' Scale:  ',,90.0,7) 

call  STRING_LENGTH (scale_char, slen) 

call  symbol (0.50,10.0, .12, scaie_char , , 90 . 0 , slen) 

call  symbol ( 0 . 75 , 9 . 25 , . 12 , ' Zone :  ',,90.0,6) 

2one_char  =  zone_table ( zone) 

call  STRING_LENGTH ( 2one_char , slen) 

call  symbol (0.75,10.0, . 12 , zone_char, ,90.0, slen) 

write  (srow_char , ' (i3 ) ' )  srow 
write  (nrow_char, ' (i3) ' )  nrows 

stringl  =  'Start  row:  ' //srow_char/ / '  Number  of  rows  ' 

//'plotted:  '//nrow_char 

call  STRING_LENGTH( stringl, slen) 
call  symbol (1.0,9.25, .12, stringl , , 90 . 0 , slen) 

write  {scol_char, ' (i3) ' )  scol 
write  (ncol_char, ' ( i3 ) ' )  ncols 

stringl  =:  'Start  col:  ' //scol_char// '  Number  of  cols  ' 

//'plotted:  '//ncol_char 

call  STRING_LENGTH( stringl, slen) 
call  symbol (1.25,9.25, .12, stringl , , 90 . 0 , slen) 


return 

end 


